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We implemented a variation after projection (VAP) algorithm based on a triaxially deformed 
Hartree-Fock-Bogoliubov vacuum state. This is the first projected mean field study that includes 
all the quantum numbers (except parity), i.e., spin (J), isospin (T) and mass number {A). Systematic 
VAP calculations with JTA-projection have been performed for the even-even sd-shell nuclei with 
the USDB Hamiltonian. All the VAP ground state energies are within 500 keV above the exact 
shell model values. Our VAP calculations show that the spin projection has two important effects: 

(1) the spin projection is crucial in achieving good approximation of the full shell model calculation. 

(2) the intrinsic shapes of the VAP wavefunctions with spin projection are always triaxial, while 
the Hartree-Fock-Bogoliubov methods likely provide axial intrinsic shapes. Finally, our analysis 
suggests that one may not be possible to associate an intrinsic shape to an exact shell model wave 
function. 

PACS numbers: 21.60. Jz,21.60.Cs,21.10.Hw 


I. INTRODUCTION 

Hartree-Fock-Bogoliubov (HFB) method has been very 
successful in describing the global properties of the 
ground states throughout the whole nuclear region. As 
a mean field method, HFB breaks the symmetries of the 
nuclear system, and can be used to study the intrinsic 
shapes. The HFB calculations with Gogny force show 
that almost all the calculated 1712 nuclei have axially 
symmetric HFB minima[l|. 

Projection can be done on a HFB vacuum to recover 
the symmetries that the Hamiltonian obeys. To test the 
quality of the projected wavefunctions, one can compare 
them with the exact shell model ones using a common 
Hamiltonian. HFB and variation after projected HFB 
calculations with shell model Hamiltonians have been re¬ 
ported by several authors For those calculations 

without projection, the HFB vacuum states are often as¬ 
sumed to be axially symmetric Q. Indeed, we will see 
below that all the calculated HFB minima in sd-shell 
nuclei, except ^“^Mg, are exactly axial with the USDB 
Hamiltonian Q. 

However, if one performs the variation of the projected 
HFB vacuum, usually called variation after projection 
(VAP) 0, it is likely that the intrinsic shape may changes 
due to the inclusion of beyond mean field correlations. 
One typical example is the ground state of ^^Mg,^hich 
is predicted to be spherical at the mean field level Q , but 
it turns out to have a quadrupole deformation when the 
correlations associated with the restoration of the broken 
rotational symmetry are considered Another exam¬ 
ple is ®®Ni, whose ground state is spherical at the mean 
field level, but is slightly deformed when performing the 
projected energy surface calculation M- 

Moreover, the triaxial (7) degree of freedom plays im¬ 
portant roles on the low-lying collective dynamics in this 


mass region [TJ. In ^'^Mg the possibility of the triaxial 
deformation in the ground state was discussed for decades 
[T^ - [T^ . and it is still being used as the testing ground 
for modern theories involving angular momentum (spin) 
projection [IM3- 

In this work, we perform VAP calculations of the even- 
even sd-shell nuclei using the USDB Hamiltonian. Here, 
we allow the 7 degree of freedom in the HFB transfor¬ 
mation. The shell model Hamiltonian conserves the spin 
(J), isospin (T), as well as the mass number (A). Hence 
a complete projection should recover all J, T, and A 
quantum numbers. This is generally very much time- 
consuming because of the 7-dimensional integration ( 3 
for J, 3 for T, and 1 for A). Presently, we can only carry 
out such extensive studies in the sd shell. For efficiency, 
we use the new techniques of Refs. [Tsl - i^ to evaluate 
the kernels for projections. 

II. THE VAP METHOD 

From a randomly chosen HFB vacuum state |$o), one 
can construct a new HFB vacuum state |$) using the 
Thouless theorem Q. Namely, 

1 $) ( 1 ) 
where d is a skew symmetric matrix, and M is the normal¬ 
ization factor. The triaxiality of the HFB vacuum can be 
treated similar to Ref. [2l| so that the Q2±i components 
of the quadrupole moment vanish. 

Projecting |$) onto good quantum numbers J, T, and 
A one gets the so called JTA-projection (similarly, TA- 
projection for T,A, etc.). One can evaluate the JTA- 
projected energy, 

EjtA = X fMMTfKKT{^\HPMKPMTKTP^\^)i‘^) 

MKMtKt 
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where ^MtKt spin, isospin and 

mass number projection operators, respectively. The 
isospin projection operator is similar to the spin projec¬ 
tion operator but in the isospin space. Ejta and the 
corresponding coefficients JkKt are obtained by solving 

^ ($1(77 - EjTA)PliKPllrKrP^\^)fKK^ = 0, (3) 

KKt 


with /kKt satisfying 

fMMTfKKT{^\PMKPMTKTP^\^) ~ W 

MKMtKt 


One can also calculate Eta with T^l-projection by sim¬ 
ply removing the spin projection from Eqs. (HI. For the 
^-projection, the corresponding energy, Ea, is reduced 
to 


^ 

Without any projection, we define 

Ehfb = ( 6 ) 

It is natural that one may consider the neutron (N) 
and proton {Z) projection, as has been done in Refs, 
id. However, this is essentially the same as the TzA- 
projection (Tz = {N — Z)j2). Here, we prefer to take 
TH-projection to recover the total isospin symmetry. 

VAP calculations can be performed by changing the d 
matrix in Eq. ©. Here, we impose the following restric¬ 
tions for the d matrix: ( 1 ) d is real, ( 2 ) keeping the time 
reversal symmetry, and (3) no mixing between neutron 
and proton in the HFB transformation. Therefore the 
total number of free VAP parameters for sd-shell is re¬ 
duced to Nvap =42. In practice we start with d = 0 
and with Nilsson-|-BCS vacuum states |<I>o) obtained with 
randomly chosen quadrupole parameters [T(TI |. 

To extract the intrinsic shape, the quadrupole moment 
and the triaxial degree of freedom, Q and 7 , are defined 
such that 



Tl\ 

/ IOtt 

(7) 

(5 cos 7 = 



{'t|\ 

/ IOtt 1 , 

( 8 ) 

Q sin 7 = 



where b is the harmonic oscillator length. |'I') refers to 
an intrinsic state, which may have different forms. Ex¬ 
plicitly, we define, 

(1) Qhfb and '^hfb for |'I') = |$), 




= , and 


( 2 ) Qa and 7 ^ for |vl-) = 

(3) Qta and jta for |4') = Y,Kt 


TABLE I: Results of the VAP-A calculations for We 

perform the VAP calculations for several times. Each time we 
start with different |$o) states. The numbers in the first col¬ 
umn denote different |'l>o) states. The second column shows 
the converged energy E^. Quantities in other columns are cal¬ 
culated with the converged |<I>) vacua. Energies are in MeV. 


|4>o) 

E\ 

Qa 


Ehfb 

Qhfb 

7nFs(°) 

(A) 

1 

-81.358 

18.284 

10.05 

-81.008 

18.005 

9.46 

8.110 

2 

-81.358 

18.284 

130.05 

-90.178 

18.371 

128.94 

9.013 

3 

-81.358 

18.284 

-109.95 

-82.684 

18.120 

-110.61 

8.259 

4 

-81.358 

18.284 

10.05 

-79.720 

17.905 

9.05 

8.000 


III. VAP CALCULATIONS FOR 

When performing the energy variation, one may find 
that there might be more than one energy minima. 
Therefore, the energy variation should be calculated sev¬ 
eral times with different starting |$o) states which are 
randomly chosen. We then identify the lowest minimum, 
and denote it with E* . Here and below, we only discuss 
the results corresponding to E*. 

In the present work, we adopt the USDB Hamiltonian 
0 . The HFB energy for ^^Mg is E*hpb = -80.965 MeV 
with the constraints ($|iV|$) = N and ($|Z|<I>) = Z. 
This is the only sd-shell nucleus for which the HFB cal¬ 
culation gives a non-axial shape with = 18.659 

and ')hpb = 11.96° (here and below the Q* and 7* 
are the shape parameters that can be associated with 
the absolute minimum for some VAP choice). Let’s 
first do the simplest VAP with only A-projection (called 
VAP-A). Since the particle number is already projected 
out, it might be unnecessary to impose a constraint to 
the average particle number of the HFB vacuum. To 
check this conjecture, we start from several different |$o) 
states and perform VAP-A. The results for few selected 
|$o) choices are shown in Table ID One can see that 
the VAP-A energies are identical {E*^ = —81.358 MeV). 
However, the corresponding Ehfb, Qhfb, Ihfb and 
{A) = ($|A|<I>) appear randomly. This means that the 
converged vacua |d)) are not unique, but correspond to 
the same A-projected state. The Qa values are the same, 
and although the 7 a values look different, they actually 
represent the same shape but with different orientations. 
Therefore, one can adopt the values Q\ = 18.284 and 
7 ^ = 10.05 to define the shape of the VAP-A mini¬ 
mum. If one imposes (A) = A = 8 , we still have 
E\ = —81.358 MeV, now the converged |4)) vacuum be¬ 
comes unique, with Ehfb = —79.720, Qhfb = 17.905, 
and ■yHFB = 9.05° (see the last line in Table|T|. However, 
for the VAP with TA-projection, the situation becomes 
a little different. 

VAP calculations with TA-projection (called VAP-TA) 
are listed in Table UTl Unlike VAP-A, even if one imposes 
(A) = A = 8 for ^^Mg, the converged |$) is still not 
unique as the Ehfb energy appears randomly. Moreover, 
the Ea energy is not unique either. Interestingly, after 
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TA-projection, those different |$) vacuum states have ex¬ 
actly the same projected energy = —82.831(MeV) 
and the same Qta = 17.295. Therefore, we can asso¬ 
ciate the shape parameter corresponding = Qta = 
17.295 to this projected minimum. Similarly, we found 
(after rotation) 7 ^^ = 7 ta = 0.09°, which describes an 
almost axial-shape. One can conclude that only Qta 
and 7 ta are meaningful in describing the shape of the 
VAP-TA projected state. 

A complete symmetry restoration is the JTA- 
projection. VAP results with JTA-projection (called as 
VAP-JTA) are shown in Table IIIII All the converged 
energies are —86.919 MeV, significantly closer to 
the shell model result Esm = —87.105 MeV. Again, both 
Ea and Ehfb in Table m can not be uniquely deter¬ 
mined, even if one enforces the (A) = A constraint. For¬ 
tunately, with the additional spin projection, all Eta 
values are found to be —79.879 MeV, and similarly the 
corresponding shape is described by Qta = 19.057 and 
7ta = 16.96°. Therefore, the quantities that can be asso¬ 
ciated with the shape of VAP-JTA wavefunction should 
also be Q*jta — Qta = 19.057 and ')*jta — Ita = 
16.96°. 

One can study the shape evolution of ^^Mg from HFB 
to VAP-JTA. In VAP-TA, Qta looks smaller than Q^fb 
in HFB, and ■Jta tends to be close to zero (axial shape). 
However in VAP-JTA, Q*jta 1® larger than the Q*jjpb la 
HFB, and 'y*jTA tends to describe a triaxial shape. This 
triaxiality in VAP-JTA, in comparison with VAP-TA, is 
likely caused by the spin projection. To determine if this 
phenomenon is more general, we performed systematic 
VAP calculations for a larger number of even-even sd- 
shell nuclei. 


IV. VAP CALCULATIONS FOR EVEN-EVEN 
sd-SHELL NUCLEI 

VAP calculations have been performed for the ground 
states of even-even sd-shell nuclei. The calculated ener¬ 
gies relative to the shell model ones are shown in Figure 
[T^. The numerical results are given in Table IIVI Here, 
we didn’t include the Oxygen isotopes and the A = 20 
isotones because their VAP-JTA energies are exactly the 
same as the shell model results {Esm)- This special case 
is discussed below. The VAP-JTA energies are much 
lower than those of HFB and VAP-TA. Moreover, The 
VAP-JTA energies for ^°Ne, ^®Ne, and Ar nuclei are ex¬ 
actly the same as the shell model results (see also Figure 
This can be understood by comparing the number 
of VAP parameters, Nvap, with the shell model dimen¬ 
sion, Njt (the total number of the independent basis 
states with good JT). Here, Nvap = 42. The Njt val¬ 
ues with J = 0 and T = 0 for both ^°Ne and ^®Ar are 
only 21. For ^®Ne, Njt for J = 0 and T = 4 is 43. It 
looks that when Njt is less than, or close to Nvap, then 
the VAP-JTA energy is likely to be the same as the shell 
model one. Indeed, for all even-even oxygen isotopes and 
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FIG. 1: (Color online) (a) Calculated ground state energies 
relative to the shell model results, Esm- (b) Relative VAP- 
JTA energy, E*jta ~ Esm, versus the shell model dimension, 
Njt, in JT subspace. 




FIG. 2: (Color online) Quadrupole moments Q and 7 values 
for (a) HFB vacuum states, (b) VAP-TA states, (c) VAP- 
JTA states and (d) VAP-HF states which is based on a slater 
determinant. 


for the N = 20 isotones, for which Njt < Nvap, we 
have obtained E*jrp^ = Esm- In Figure [T) d, one can also 
see that the energy difference — Esm increases with 
Njt- The largest — Esm = 0.446 MeV is obtained 
for ^®Mg, corresponding to the largest Njt = 1132. 

The quadrupole moment and the 7 degree of freedom 
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TABLE II: Similar to Table. |T]but for the VAP-TA calculations. (A) = 8 is imposed. 


|3>o) 

AJ.^(MeV) 

Qta 

7ta(°) 

AA(MeV) 

Qa 

7a(°) 

A frFs(MeV) 

Qhfb 

7irFs(°) 

(A) 

1 

-82.831 

17.295 

-119.91 

-75.826 

16.376 

-118.62 

-74.921 

15.755 

-118.23 

8.000 

2 

-82.831 

17.295 

0.09 

-74.402 

16.167 

2.47 

-73.909 

15.563 

3.06 

8.000 

3 

-82.831 

17.295 

120.09 

-76.633 

16.526 

120.09 

-75.525 

15.897 

120.08 

8.000 


TABLE III: Similar to Table. |T]but for the VAP-JTA calculations. (A) = 8 is imposed. 


|3>o) 

AjFA(MeV) 

ATA(MeV) 

Qta 

7ta(°) 

AA(MeV) 

Qa 

7a(°) 

E ffFs(MeV) 

Qhfb 

'yHFBO 

(A) 

1 

-86.919 

-79.879 

19.057 

-16.964 

-75.600 

17.482 

-20.225 

-73.781 

16.230 

-23.772 

8.000 

2 

-86.919 

-79.879 

19.057 

-16.963 

-75.641 

17.510 

-20.119 

-73.830 

16.264 

-23.604 

8.000 

3 

-86.919 

-79.879 

19.057 

-16.963 

-75.644 

17.520 

-20.068 

-73.845 

16.281 

-23.506 

8.000 


can be extracted using Eqs. dZ]) and (HI). In Fig. [2^, 
the values in HFB are either 0° or 60°, except 

')hfb — 12 ° for ^^Mg, thus supporting the conclusion 
that HFB likely presents axially deformed shapes. In 
Fig. [2 )d, the shapes in VAP-TA calculations still remain 
axially symmetric, except for ^®Mg, which has 7 ^^ = 
25.7°. Quite differently, the values in the VAP-JTA 
calculations (Fig. show that all these nuclei are non- 
axial without exception. Comparing these results with 
those of Fig. [5^, one can conclude that the triaxiality in 
VAP-JTA is definitely a beyond mean-field effect, which 
is likely to be a universal phenomenon. Fig. however, 
excludes the possibility that the isospin projection and 
the mass projection lead to triaxiality. Thus, the only 
possible cause of the triaxiality is the beyond mean-field 
spin projection. 

To study directly the effect of spin projection, one can 
start from a Hartree-Fock (HF) Slater determinant (SD) 
and perform VAP calculations with only spin projection 
(called VAP-HF). The converged energies, E*ppp, rela¬ 
tive to Esm, are shown in Fig. [1^. The results show 
that VAP-HF is better than VAP-TA, and quite close 
to the VAP-JTA. The quadrupole moment Q*ppp and 
'y’pHF corresponding to Eppp can be calculated using 
Fqs. d?]) and dH]) with Idt) replaced by the converged SD. 
These quantities are uniquely determined, and are shown 
in Fig. [2ji. Again, all the Jppp values are distributed 
in the interval (0°, 60°), which is very similar to Fig. [2t. 
Therefore, we could conclude that VAP results that in¬ 
clude spin projection can always be associated with in¬ 
trinsic states having triaxial deformation. 

One more interesting phenomenon, however, is re¬ 
lated to the VAP-JTA calculations for ^°Ne, ^®Ne, and 
®®Ar. We have shown above that the energies of 

these nuclei are the same as the exact shell model re¬ 
sults. Surprisingly, the corresponding Qta and 7 ta val¬ 
ues are not unique, which is quite different from other 
nuclei with Ejp^^ > Esm- For example, the results for 
20 Ne are shown in Table El With the same converged 
Ejta = —40.472MeV, one can clearly see that starting 
with different initial states |$o), the result for Qta and 
7 ta could be different. These results indicate that it may 


not be possible to associate an unique intrinsic deforma¬ 
tion with an exact eigenstate of the Hamiltonian. 


V. SUMMARY 

We implemented an algorithm that performs variation 
after projection (VAP) on spin, isospin, and mass number 
of a triaxially deformed Hartree-Fock-Bogoliubov vac¬ 
uum state. This is the first projected mean field study 
that includes all these quantum numbers. 

We start from a randomly chosen HFB vacuum state 
and carry out VAP calculations for ^^Mg in sd-shell with 
various projections. In the VAP-A case the converged 
solution is independent of the Fermi level (chemical po¬ 
tential) . Although the associated HFB vacuum does not 
have definite quadrupole moment Qhfb and triaxial de¬ 
formation parameter jhfb, one can use the unique Qa 
and jA to describe the intrinsic deformation of the VAP- 
A state. Similarly, in the VAP-TA calculations, Qa and 
jA can not be uniquely determined, but Qta and 7 ta 
are unique and can be associated with the intrinsic defor¬ 
mation of the VAP-TA state. It is not possible to directly 
define deformation parameters Q and 7 for the VAP-JTA 
wave function, which has the symmetries fully restored, 
but the Qta and 7 ta calculated with the VAP-JTA vac¬ 
uum state 1 4)) are also unique, and can be associated with 
the intrinsic deformation of the VAP-JTA state. 

Systematical VAP calculations of even-even sd-shell 
nuclei have been performed using the USDB Hamilto¬ 
nian. The VAP-JTA energies, E*jp^, are very close to 
the shell model results, Esm- Moreover, the relative en¬ 
ergy, Ejp^ — Esm, increases with the shell model dimen¬ 
sion Njt- The shapes described by the HFB minima are 
always axial. However, with spin projection VAP calcu¬ 
lations always produce triaxial shapes. We believe that 
such triaxiality is an universal phenomenon caused by 
the beyond mean-field dynamic correlations. Finally, we 
show that those VAP-JTA states reaching the exact shell 
model results do not have clearly defined intrinsic shapes. 
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TABLE IV: Converged energies and associated shape parameters for even-even sd-shell nuclei calculated with the USDB 
Hamiltonian. 


Nucleus 

1 VAP-JTA 

1 VAP-TA 

1 HFB 

1 VAP-HF 


Njt 

Esm 

IP* 

^JTA 

Q JTA 

Ijta 

Eta 

Qta 

Ita 

Ehfb 

Qhfb 

Ihfb 

TP* 

^PHF 

Qphf 

Iphf 

"‘’Ne 

21 

-40.472 

-40.472 

- 

- 

-37.069 

14.7 

0.0 

-36.404 

15.3 

0.0 

-40.265 

13.861 

3.551 

^^Ne 

148 

-57.578 

-57.501 

12.1 

13.8 

-54.572 

15.8 

0.0 

-53.474 

16.5 

0.0 

-56.958 

15.675 

8.632 

^^Ne 

287 

-71.725 

-71.570 

11.0 

30.1 

-68.084 

10.1 

60.0 

-66.402 

12.0 

0.0 

-71.037 

13.449 

32.786 

2 ®Ne 

191 

-81.564 

-81.465 

9.2 

28.4 

-78.949 

8.6 

0.0 

-77.518 

8.3 

0.0 

-80.988 

9.760 

17.265 

2 ®Ne 

43 

-86.543 

-86.543 

- 

- 

-84.920 

7.0 

60.0 

-83.949 

7.1 

0.0 

-86.294 

9.848 

23.934 

'^Mg 

325 

-87.105 

-86.919 

19.1 

17.0 

-82.831 

17.3 

0.0 

-80.965 

18.7 

12.0 

-86.636 

19.165 

16.427 

2® Mg 

1132 

-105.521 

-105.075 

15.8 

28.7 

-100.648 

13.8 

25.7 

-98.992 

15.9 

60.0 

-104.264 

16.238 

32.331 

Mg 

874 

-120.500 

-120.205 

14.5 

20.2 

-117.091 

14.4 

0.0 

-115.625 

15.1 

0.0 

-119.354 

16.306 

19.835 

®°Mg 

191 

-130.474 

-130.400 

10.4 

20.0 

-128.035 

10.3 

0.0 

-126.735 

10.9 

0.0 

-129.926 

11.864 

27.322 

"»Si 

839 

-135.860 

-135.539 

16.1 

58.6 

-131.501 

17.8 

60.0 

-130.021 

19.8 

60.0 

-134.617 

17.116 

58.038 

®°Si 

1132 

-154.754 

-154.402 

14.3 

47.0 

-150.380 

10.6 

60.0 

-148.475 

14.5 

60.0 

-153.777 

14.633 

46.167 

®2Si 

287 

-170.519 

-170.373 

12.5 

58.2 

-167.721 

10.6 

60.0 

-166.344 

12.4 

60.0 

-169.996 

12.175 

52.023 

32s 

325 

-182.452 

-182.234 

15.1 

33.3 

-179.925 

0.6 

60.0 

-176.393 

0.0 

0.0 

-181.856 

14.681 

33.188 

34s 

148 

-202.504 

-202.380 

10.6 

53.0 

-200.331 

0.0 

0.0 

-198.493 

0.0 

0.0 

-202.039 

11.496 

50.992 

®6Ar 

21 

-230.277 

-230.277 

- 

- 

-228.355 

0.0 

0.0 

-226.611 

13.2 

60.0 

-230.112 

12.068 

52.293 


TABLE V: VAP results with JTA projection for ^°Ne. 


l-J’o) 

E}^^(MeV) 

ATA(MeV) 

Qta 

JtaC) 

1 

-40.472 

-28.284 

6.314 

-45.134 

2 

-40.472 

-30.468 

11.873 

-124.746 

3 

-40.472 

-27.932 

9.876 

2.592 
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